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Supersymmetric S'O(IO) Grand Unified Tlieories with Yukawa unification represent 
an appealing possibility for physics beyond the Standard Model. However Yukawa uni- 
Q>^ ' fication is made difficult by large threshold corrections to the bottom mass. Generally 

, one is led to consider models where the sfermion masses are large in order to suppress 

' these corrections. Here we present another possibility, in which the top and bottom 

GUT scale Yukawa couplings are equal to a component of the charged lepton Yukawa 
matrix at the GUT scale in a basis where this matrix is not diagonal. Physically, this 
weak cigenstatc Yukawa unification scenario corresponds to the case where the charged 
leptons that are in the 16 of SO (10) containing the top and bottom quarks mix with 
their counterparts in another 50(10) multiplet. Diagonalizing the resulting Yukawa 
matrix introduces mixings in the neutrino sector. Specifically we find that for a large 
, region of parameter space with relatively light sparticles, and which has not been ruled 

^ • out by current LHC or other data, the mixing induced in the neutrino sector is such 

5^ , that sin^ 2023 ~ 1 , in agreement with data. The phcnomenological implications are 

analyzed in some detail. 



1 Introduction 



The idea of unification has been central to the development of modern physics. Supersym- 
metric (SUSY) Grand Unified Theories (GUTs) unify the gauge groups of the standard 
model as well as fermions and bosons, in addition to potentially providing explanations 
of a variety of theoretical and observational problems or phenomena [Il[2] . The choice of 
SO{10) allows for the unification of all matter superfields (of each generation) into a 16- 
dimensional spinor representation, provided we add a superfield containing a right-handed 
neutrino to complete the multiplet [SHS] . An advantage of including right-handed neutrinos 
is that it allows for the natural incorporation of the see-saw mechanism to generate neu- 
trino masses [6l410j. Finally, in the simplest SO{W) SUSY GUT models one has Yukawa 
unification, by which one means that e.g. 

yt = Vb = Vt, (1) 

where yt, Ub, and yr are the top, bottom, and tau Yukawa couplings at the GUT scale. 
(Simple Yukawa unification is not phenomenologically viable for the first and second gen- 
eration.) 

A great deal of work has gone into studying Yukawa unification (see for example |11H72| ) 
and in particular understanding whether Yukawa unification is theoretically consistent and 
compatible with observation. This is an important and non-trivial question, due in part 
to the large weak scale threshold corrections (corrections that occur when one matches the 
SM to the MSSM), especially to the bottom mass |69fl71] . These threshold corrections 
make unification of third generation Yukawa couplings at the GUT scale dependent on the 
SUSY spectrum. Thus it is only for certain values of the soft SUSY breaking parameters 
that Yukawa unification is obtained. 

To be more specific, the parameters of the SUSY 50(10) models considered here are 

7711/2, ""T-ie) "iio, M|,, tan/3, and sign(/_i). (2) 

We assume the following boundary conditions at the GUT scale, 

mQ=m\ = ml = ttiIq + Ml) (3) 

ml = ml = m?6 - 3M|, (4) 

= + (5) 

"^k,. = ^lo T 2M|,. (6) 

Here, the parameters ttiq, m\, rnfj, fn'jj, and mj^ refer to the squared SUSY breaking 
masses for the left-handed squarks, right-handed sleptons, right-handed up-type squarks, 
right-handed down-type squarks, and left-handed sleptons respectively. Of course, to ob- 
tain physical sfermion masses, one must find the eigenvalues of the relevant sfermion mixing 
matrix. Similarly is the squared mass of the right-handed sneutrino, and give 
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the squared Higgs mass terms for the Higgs doublets which couple to up-type quarks, and 
down- type quarks and leptons, respectively. The parameters mi/2, mig, and mio are the 
GUT scale masses of the gauginos, the sfermions, and the Higgs soft terms, respectively. 
D-term contributions to GUT scale Higgs and sfermion soft terms are parameterized by 
M|,. These terms are necessary to ensure radiative electroweak symmetry breaking [72| . 
While in principle, either sign of M|, is possible, we will only consider M|, > as this 
is the sign which lowers fnjj^ with respect to rn'jj^ and hence aids in allowing radiative 
electroweak symmetry breaking when m, ~ yt- Yukawa unification only occurs for special 
regions in the parameter space defined above. Generally these are regions with quite heavy 
scalars and fairly light gauginos. Typically there is also an inverted mass hierarchy, which 
is radiatively generated [75H75] . 

There are several difficulties with this parameter space. First, the relatively light 
gluino is disfavored by the non-observation of any excess in jets plus missing energy at 
the LHC |79H100] . Secondly, in this region of parameter space, the mostly bino Lightest 
Supersymmetric Particle (LSP) will tend to have too small of an annihilation cross section 
to be a (thermal) dark matter candidate |42tl45^ [50]. Finally, raising the sfermion masses 
tends to suppress contributions to the g — 2 of the muon, which therefore remain in tension 
with the experimental value [lOlj . Thus, we are interested in finding a way to preserve 
the elegance of Yukawa-unified 50(10) models while also allowing for somewhat heavier 
gauginos and lighter sfermions. 

Our approach will be to unify the top and bottom Yukawa couplings at the GUT 
scale with the third generation diagonal element of the charged lepton Yukawa matrix, but 
to introduce an arbitrary mixing angle between the second and third generation charged 
leptons |102j . Additionally, we will assume that the right-handed neutrino in the third 
generation 16 mixes with additional SU (5) singlet states at the scale at which 50(10) is 
broken to SU{5); this effectively lowers the neutrino Yukawa coupling. 

Assuming that these are the only large mixings, the charged lepton mixing angle will 
translate to the 623 of the PMNS matrix |103pi04j which governs the neutrino sector. Thus 
we can view the parameter which we have added to the standard Yukawa-unified SUSY 
picture as postdicting (or being constrained by) large mixing in the neutrino sector. This 
scenario could be embedded into some specific GUT realization. We will not consider a 
specific embedding here; a discussion of this possibility was presented in |102| . 

In Section [2] we will describe this alternative scenario, in which second and third gener- 
ation mixing in the charged lepton sector is used to accommodate third generation Yukawa 
unification, in more detail. The consequences of this mixing for the neutrino sector will 
be discussed in Section [3l In Section U] we will describe a scan of the parameter space de- 
scribed in Eq. ([2]). In particular we will note the preferred values for the mixing angle and 
its consequences for the neutrino sector in Section [5l while we will discuss the collider and 
dark matter phenomenology of parameter space points allowed by the scan in Section [6l 
We present our conclusions in Section [71 
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2 The Scenario 



In the absence of threshold corrections, Yukawa unification is relatively easy to achieve. In 
fact values of the top mass close to that which was in fact measured had been predicted 
in Yukawa unified scenarios (see e.g. [T11[T21[SU1[7T]). Due to the threshold corrections 
mentioned above, however, the bottom Yukawa at the weak scale is different from what it 
would be naively. These corrections lower the bottom mass when //M3 > 0, as is suggested 
by the observed discrepancies in 5 — 2 |101j . In this case, a larger value of tan/3 than the 
~ 55 expected from taking the ratio of the running top and bottom masses at the weak 
scale is required for the unification of yt and yb at the GUT scale. The corresponding 
threshold corrections to the tau mass tend to be smaller than those to the bottom mass, 
due to the absence of SUSY QCD contributions. As a result, for the values of tan /? required 
for top-bottom unification with threshold corrections, the tau Yukawa at the GUT scale 
tends to be larger than the (unified) top and bottom Yukawa couplings. 

Here we consider a scenario where instead of unifying the tau Yukawa with the the top 
and bottom Yukawas at the GUT scale, we take the second and third generations of the 
charged lepton Yukawa matrix to be described by 

UMovt)=(''' (7) 

V ^ yt/bj 

where y^/^ is the (unified) top and bottom Yukawa coupling at the GUT scale. We have 
limited ourselves to a symmetric matrix for simplicity. In this paper, we will consider only 
the second and third generations. We do not expect that considering all three generations 
would qualitatively affect our results, though future research will evaluate this claim. 

The form of the Yukawa matrix in Eq. d?]) describes a scenario where the third genera- 
tion Yukawa coupling at the GUT scale is yt/i, for an entire 16 of S'O(IO), but the charged 
leptons in this 16 mix with the charged leptons in the second generation 16. To be precise, 
both the left and right-handed charged leptons in the "third generation" 16 mix with their 
counterparts in the second generation 16 with the same mixing angle Q. 

As a result of the assumption that the mixing angle is the same for both charged lepton 
states, the transformation of the Yukawa matrix from the basis in which it is diagonal is 
simply that resulting from a rotation. Thus, due to the invariance of the trace and the 
determinant under rotations, the parameters in Eq. ([7]) are given by 

y22 = yr + ytM- Vt/b (8) 

and 

= {Vt/b)y22 - VrVfi- (9) 

In these equations, as in Eq. ([7]) above, all Yukawas are to be evaluated at the GUT scale. 
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We can also write the non-diagonal charged lepton Yukawa matrix of Eq. ([7]) in terms 
of Ut, y^, and the mixing angle 0, as follows: 



Vt + Vim fl 0\ Vt- Vfi f- cos 26 sin2e\ , , 

2 VO V 2 V sin 26 cos26y' ^ ' 

where 6 describes the mixing of the charged leptons in the 16 with the top and bottom 
with the charged leptons in another multiplet. We see immediately that 

where we consider the particular quantity sin^ 26 here due to its significance in the induced 
mixings in the neutrino sector, which are discussed in Section [3] below. 
Using Eq. ^ and Eq. ([9]), Eq. (fTT]l can be written explicitly as 

. 2 _ ^ivt/biyr + Vt,- Vt/b) - Vry^i) 

sin Z\y — . \^^) 

[Vt - VfiY 

We now define 

R = ^ (13) 
Vt/b 

for convenience. We note that in the case where the top and bottom Yukawas are equal at 
the GUT scale, this corresponds to the R in e.g. [54] when yr > yi/i,, as is in fact obtained 
in this scenario. Ignoring y^, as compared with yr, Eq. (fT2]) can be re-expressed as 

sm^2Q^4^^ (14) 

Thus, we will find sin^ 26 ~ 1 when i? ~ 2. In fact, this result is quite stable for i? ~ 2; to 
be more precise, if -R = 2 + e, then Eq. [TH becomes 

sin^ 26 « -^P^ ^ 1 - 6^4. (15) 

Therefore, if we write our charged Yukawa matrix at the GUT scale in the form given in 
Eq. ([7]) or Eq. (jlOp . where the Yukawa matrix has a diagonal element equal to the top and 
bottom Yukawas at the GUT scale, then the mixing will be maximal when the tau Yukawa 
at the GUT scale is approximately twice the (unified) top and bottom Yukawa coupling. 

The value of 6 obtained above will, in general, evolve as the charged lepton Yukawa 
matrix is run from the GUT or Majorana mass scale to the weak scale. However, in 
the limit where neutrino Yukawa couplings vanish, every term in the one or two loop /3 
functions |105| for Ye is either proportional to Yg, YgYlYg, or YeYlYeY^Ye. Writing 
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Ye = C/^^DgC/, we find that we can write YeY^Ye = C/^DeD*DeC/, etc. Thus we can write 
the Renormahzation Group Equations (RGE) in the schematic form: 

'^^ = u^mM (16) 

where f(De) is the right hand side of the RGE in terms of diagonal Yukawa matrices. 
Since the RGE (obviously) take the same form for diagonal Yukawa matrices as for general 
Yukawa matrices, we have that 

U^^U = U^fiT>,)U, (17) 

and, result, 

dW . U 

—BeU + U^'De— = 0. 18) 
dt dt ^ ^ 

Plugging in the general form of the rotation matrix for U and our specific values for the 

entries of De, we obtain 



fde\ /sin 26 cos 26 \ ^ 

[^)(yr-y.)[,,,2e -sin2ej=°' ^'^^ 



which, given the non-degeneracy of i/r and y^i, can only hold if 

Thus, in the limit where neutrino Yukawa couplings are zero, the value of Q (or functions 
of 0) obtained, e.g. in Eq. (fH|) is the value at any scale, not just at the GUT scale. While 
neutrino Yukawa couplings are non-zero, and in fact generally large in Yukawa unified 
models (at least for the third generation), as mentioned before we will assume a dilution 
of neutrino Yukawa couplings via the mixing of the right-handed neutrino in the third 
generation 16 with other states. Therefore the RGE effects on the mixing angle G become 
negligible. 



3 Neutrino Sector 
3.1 Mixing 

As noted above, in this scenario, the third generation fermions are part of a 16 of 50(10) 
for which Yukawa unification at the GUT scale holds, but where the right- and left-handed 
charged leptons in this multiplet mix with the right- and left-handed charged leptons in 
the second generation 16. We assume that there is no similar inter-generation mixing in 
the neutrino sector. 
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As the W will only couple states within these multiplets, the left-handed neutrino in the 
16 will couple to a linear combination of the physical muon and tau. Considering instead 
the combinations of left-handed neutrino eigenstates from the two multiplets which couple 
to a given charged lepton mass eigenstate, we find 

i^fi = cos Q 1^2 — sin Q u^, (21) 

and 

i^r = sin G 1^2 + cos Q z/3 , (22) 

where we have defined the non-diagonal charged lepton matrix in terms of the mixing angle 
0, as in Eq. (jlOp above. This gives the 2 by 2 PMNS matrix U defined by 

/^M_/^cosG -sine\ /z^2\ /^,o^ 
\ur) ~ {sine cose J [usj ^ 

From this we find that 

023 = -e (24) 

Thus, in particular, we find 

sin^ 26*23 = sin^ 26. (25) 

3.2 Yukawa Couplings 

It has been noted, for example in |106fllll] . that large values of neutrino Yukawa couplings 
generically lead to large lepton fiavor violation (though for a somewhat contrary perspective 
see |112j ). This flavor violation arises from the contribution of this Yukawa term to the 
RGE running of the slepton mass matrix. As noted above, we assume a dilution of neutrino 
Yukawa couplings, e.g. through additional mixing between the right-handed neutrinos and 
unspecified additional states at the GUT scale. We demand sufficient dilution to reduce 
BR{fj, 67) to below the experimental bounds, in particular the 90% C.L. upper limit on 
this quantity, 2.4 x 10"^^ |113j ). We will examine this mixing and the resultant dilution of 
the neutrino Yukawa couplings in greater detail in future work. 

4 Scan of Parameter Space 

As a step toward determining which regions in the parameter space of this scenario are 
allowed by existing constraints, we performed a scan of the parameter space using the SUSY 
spectrum generation and analysis code SOFTSUSY 3.1.7 |114j , as well as the dark matter 
calculation code DarkSUSY 5.0.5, |115] and the code for NLO calculation of BR{b — > 57) 
in the MSSM, SusyBSG 1.4 |116j . We modified the GUT-scale boundary conditions in 
SOFTSUSY to be those of this scenario (Eq. (g]) - ([6])). 
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Mill 


Max 


^'I'PTi '\XT^ H "1" n 


mi/2 


400 GeV 


2000 GeV 


25 GeV 


miQ 


300 GeV 


2000 GeV 


20 GeV 




1.1 


1.4 


0.05 




0.275 


0.4 


0.025 




-2000 GeV 


-200 GeV 


30 GeV 


tan /3 


58 


64 


0.3 



Table 1: Ranges of parameters used in the parameter scan described in Section 4. 



In this scan, the supersymmetric parameter ^ was always chosen to be positive, in order 
to reduce the tension with the observed value of muon g — 2 [lOlj . The ranges for the other 
parameter regions in which the points are found are given in Table HI points were scanned 
using a grid search. We found that with the exception of the upper bound for mig, the 
specific values of the bounds used in the grid search were not particularly important, as 
very few points in the parameter space that satisfied the constraints we impose (detailed 
below) were near the boundary of the parameter space scanned. The upper bound on 
mie was used to restrict our focus to regions of parameter space with somewhat lighter 
sfermions. 

A given point was only considered for further study if the following criteria were satisfied 

1. Consistency of spectrum. Specifically none of the "problem flags" listed in Ap- 
pendix B of Ref. |114| are triggered. These flags are triggered if the spectrum 
calculation does not converge to the specified tolerance, if couplings become non- 
perturbative, if a Landau pole is encountered in the RGE running, if electroweak 
symmetry is not broken, if there are tachyons in the model, if the p parameter can- 
not be calculated, if the Higgs potential is unbounded from below, or if the GUT 
scale is less than 10^ or greater than 5 x 10^^ GeV. 

2. Top/ bottom Unification. The top and bottom have the same Yukawa coupling 
at the GUT scale. Specifically we demand 

I ~ I < 0.02 (26) 

yt + Vb 

where yt and yt are the top and bottom Yukawa couplings at the GUT scale. 

3. R in desired range. In the scenario described, if i? > 2, then y22 > yss- {R is 
defined in Eq. ([131) •) We have restricted our analysis to the case in which 7/33 > 7/22 
and hence R < 2. Very few parameter points were eliminated by this constraint. 

4. SUSY spectrum constraints. Charged sparticles must be heavy enough to have 
evaded detection at LEP; specifically we demand that the chargino mass be greater 
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than 94 GeV. Correspondingly, the neutrahno mass is greater than 46 GeV. In addi- 
tion, we require that the hght SM-Uke Higgs mass is greater than 114 GeV. We note 
that some of these constraints on the SUSY spectrum are conservative; our main 
concern is with determining the general properties of the allowed parameter space, 
rather than exploring ways in which naive bounds could be circumvented. 

5. Dark Matter Constraints. We demand that the LSP be a neutralino, and that 
its thermal relic density satisfy 0.1 < ^x^^ ^ 0.12, following |117j . We implement 
the bound on the spin-independent dark matter direct detection cross section from 
XENONIOO |118| . which is ~ 10~^^ cm^ for most LSP masses. Likewise, to avoid 
bounds e.g. from COUPP |119j . we demand that the spin-dependent WIMP-nucleon 
cross section be less than 10~^^ cm^. 

6. B Physics Constraints. We demand that the calculated value of BR{h — )■ 57) be in 
a range consistent with the world-average experimental value 3.55 it 0.24 it 0.09 |120j . 
Specifically, we demand that the value of BR{b — t- 57), calculated with SusyBSG 
1.4 [nS], be in the range 3.04 x 10""^ < BR{h 57) < 4.06 x 10"^. Likewise, 
we demand BR{Bs — )• fi^ < 1.2 x 10~^, for consistency with recent LHCb and 
CMS searches |121H125] for this rare decay. We do not add a separate contraint 
from BR{B — >■ tv), as the value of this quantity in the parameter space points 
which satisfy the preceding constraints, epsecially the limit on BR[Bs — )• fj,~^fi~), is 
always consistent with experiment as will be discussed in more detail below. Both 
BR{Bs I^^IJ-') and BR{B tu) were calculated with Superlso 3.2 [T26l[T27] : 
the calculations for all Higgs masses and B physics observables were verified using 
CPsuperH2.0 [128] . 



5 Values of sin^ 2^23 

The distribution of sin^ 2^23 for models obtained in the scan described above is presented in 
Fig. [TJ Note that all of the models found in the scan are consistent with the experimental 
observation that sin^ 2^23 ~ 1- We now examine why this relation holds so generally. 

The value of sin^ 2^23 is found from Eq. (jl2p . As we see from the approximation in 
Eq. (jl4p . this value is mainly a function of the (unified) top and bottom Yukawa coupling 
y^/f, as well as the tau Yukawa coupling ?/,-, defined at the GUT scale. The distribution of 
ytjii and yr are shown in Figs. [2] and [3] respectively. We note that the unified top/bottom 
Yukawa coupling is always in the range 0.55 — 0.59, while the tau Yukawa coupling ranges 
from ~ 0.9 to ~ 1.1. We note, following Eq. (|14p . that the lower end of this range 
corresponds to values of sin^ 2^23 of ~ 0.93, while the upper end of the range corresponds 
to values of ~ 1. Thus having sin^ 2^23 in the desired range is a result of the relatively 
large tau Yukawa coupling. 

We can understand the large Yukawa coupling as follows. In the absence of SUSY 
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Figure 1: Distribution of sin^ ^23 (obtained from taking 623 = —Q, and using Eq. ()12p to 
determine 0) for the models found in the scan. 
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Figure 2: Distribution of the (unified) top/bottom Yukawa couphng at the GUT scale in 
the models found from the scan. (In our scan, we demand that these top Yukawa couplings 
be unified to within 4% at the GUT scale; we take the geometric mean of these quantities 
when constructing these histograms.) 
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Tau Yukawa Couplings at GUT Scale 
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Figure 3: Distribution of the tau Yukawa coupling at the GUT scale in the models found 
from the scan. 



corrections to the bottom mass, the value of tan (3 which leads to top-bottom unification 
yields a value for the GUT scale tau Yukawa coupling which is similar to the value for 
the top and bottom Yukawa couplings. Standard Yukawa unification is not hard to obtain 
in this scenario. Adding the large negative contributions to the bottom mass raises the 
value of tan /? needed to obtain top-bottom unification. This has relatively little effect on 
the value of the unified top-bottom Yukawa coupling, but serves to raise the tau Yukawa 
coupling into the desired range, as the tau Yukawa is also enhanced by tan (3. 

Hence we would expect the (absolute value) of the fractional correction to the bottom 
mass to correlate with larger values of sin^ 2^23 ) which in fact we do find (as shown in 
Fig. HI) . Additionally, we would expect that the value of tan (3 is correlated with sin^ 2^23 • 
This expectation is confirmed in Fig. [5l Since the dominant corrections to the bottom 
mass are proportional to the gluino mass, but inversely proportional to the square of the 
sbottom mass, we expect that models with heavier gauginos, i.e. larger values of 
will have larger corrections to the bottom mass, and hence larger values of sin^2023- On 
the other hand, we expect larger values of mig, which correspond to larger values of the 
sbottom mass, to suppress the SUSY correction to the bottom mass and hence the value of 
sin^ ^23 obtained. These expectations for the dependence of sin^ 2^23 on myj and mie are 
born out in the models found by the scan, though the dependence in each case is somewhat 
weak. 



11 



sin^2e23 and the Fractional Correction to tine Bottom Mass 
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Figure 4: Distribution of threshold correction to the bottom mass as a fraction of the 
bottom mass and the value of sin^ 2^23 (obtained as above) in models found in the scan. 




Figure 5: Distribution of the values of tan/3 and sin^ 2^23 (obtained as above) in models 
found in the scan. 
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CP-Odd Higgs Mass in Models from the Scan 
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Figure 6: The distribution of the CP-odd Higgs mass is shown in this histogram. The heavy 
CP-even Higgs and the charged Higgses are nearly degenerate with the CP-odd Higgs in 
every model found by the scan, so this histogram also describes the mass distribution of 
these other Higgses. 



6 Phenomenology 

We now consider the consequences of this scenario for collider and dark matter phenomenol- 
ogy. In each case, an important difference between the scenario described here and more 
standard Yukawa-Unified SUSY scenarios, e.g. as considered in [57], is that third genera- 
tion sfermions may be relatively light in this scenario. The consequences of this difference, 
as well as the general signatures to be expected from points in the parameter space found 
by our scan (which we will refer to as "models" for convenience), will be discussed below. 



6.1 Higgs Phenomenology 

Models found by the scan all have decoupled Higgs sectors; i.e., in all models the CP- 
odd Higgs is heavy and the heavy CP-even Higgs and the charged Higgs bosons are nearly 
degenerate with it. (Specifically, 0.993 ruA < mn < ruA and 1.001 niA < 'mfj± < 1.010 uia 
in all models found by the scan.) The distribution of CP-odd Higgs masses is shown in 
Fig. [6l In every model found by the scan, the light CP-even (Standard Model-like) Higgs 
has a mass between 118 and 122 GeV; the distribution of these masses in the model set 
is shown in Fig. [71 Thus a Higgs mass between 118 and 122 GeV is a prediction of this 
scenario, up to uncertainties in the Higgs mass calculation. 
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SM-Like Higgs Mass in Models from the Scan 
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Figure 7: This histogram shows the distribution of hght, CP-even Higgs masses in models 
found by the scan. 



6.2 Flavor Physics and Precision Observables 

One can observe from Fig. [6] that the CP-odd Higgs tends to be somewhat heavy. This is 
due in part to the constraint on BR{Bs fx^ ^~), which is effectively a constraint on the 
CP-odd and charged Higgs masses. As was noted above, the heavy CP-even Higgs and the 
charged Higgs bosons are nearly degenerate with the CP-odd Higgs in all of the models 
found by the scan. The dependence of the values of BR{Bs — )• fi'^ on the CP-odd Higgs 
mass is shown in Fig. [8l While other parameters beside the CP-odd Higgs mass affect the 
value obtained for this rate, the requirement that BR{Bs fi^fj,~) < 1.2 x 10"*^ effectively 
sets a lower limit of ^ 600 GeV on the CP-odd Higgs mass. 

Fig. [9] shows the dependence of BR{B — ?• tv) on the CP-odd Higgs mass (which is 
nearly degenerate with the charged Higgs mass). The values of BR{B — ?> tv) shown in 
this figure were calculated using Superlso 3.2 |126U127] . but have been scaled so that the 
SM value is 1.22 it 0.31 x 10~^, which is calculated using the value of Vub obtained from 
inclusive decays |130j . Using the experimental average for this branching ratio from the 
HFAG cohaboration [I20l[l3l], 1.64 ± 0.34 x lO"'^, we find that the 2a lower bound on 
the ratio of the MSSM value for BR{B — >■ tv) to the SM value is 0.46. This corresponds 
to a value for BR{B — >■ tu) of 0.56 x 10~^, which is somewhat below the minimum value 
obtained for parameter space points found by our scan. Due to the l/m^± dependence of 
the supersymmetric contribution to the decay amplitude, the smallest values of this decay 
rate are obtained for the smallest values of the CP-odd mass, as it is clearly seen in Fig. [H 
In this case, BR[B — > tu) is large enough when the CP-odd Higgs mass is > 600 GeV, 
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Figure 8: This figure shows the dependence of BR{Bs — )■ ) on the charged Higgs 

mass. 

which was necessary for the BR{Bg — t- /i "*"//") constraint to be fulfihed. 

We note that there is tension between the SM theory value and the experimental average 
for BR[B — )■ Tv), and that this tension is somewhat exacerbated in the parameter space 
points found by our scan. This is in fact a generic feature of SUSY models. This tension 
is significantly increased if one uses the value of Vub obtained from exclusive decays; in 
this case the SM prediction is 0.67 it 0.15 x 10""^ |130| . and the 2a lower limit on the 
ratio of MSSM BR{B tv) to SM BR{B tv) is 0.95. This corresponds to a value of 
0.64 X 10~^. However, if 0.67 ± 0.15 x 10~^ is used for the SM theory value, the values 
displayed in Fig [9] would need to be scaled down by 1.22/0.67 = 1.8 (as it is the ratio of 
MSSM to SM rate which has been calculated). Thus this constraint would in fact rule out 
all points found by our scan, effectively demanding a CP-odd Higgs mass > 1300 GeV. 

The distribution of the branching ratio for h ^ \s shown in Fig. [TOl We note that the 
models found by the scan predominantly have values of this quantity, which like the SM 
value (3.15 =b 0.23 |129j ) are somewhat below the world-average experimental value |120j . 
though it should be noted that the theory error is comparable to this discrepancy. The 
SUSY contribution to the anomalous magnetic moment of the muon is shown in Fig. [TTl 
We note that the value of this quantity obtained in these models generally reduces the 
tension between the SM value and experiment (the experimental value is higher by (22.4 it 
10) X 10"^° to (26.1 ±9.4) X 10-^° [M] depending on the precise method used to calculate 
the SM contribution). However the size of the SUSY contribution in models found by the 
scan is insufficient to account for the experimental central value; generic models found by 
the scan have a SUSY contribution of ~ (4 — 8) x lO"^'^. An extension of scenario presented 
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Figure 9: This figure shows the dependence of BR{B — ?> tv) on the charged Higgs mass. 

here, in which the first and second generation sfermion mass terms do not have the same 
boundary conditions as the third generation sfermion mass terms, might allow for g — 2 
values closer to the observed value. 

6.3 Gluino Production and Cascade Decays 

The gluino mass and the decay modes of the gluino play a large role in determining what 
signatures a given SUSY scenario will have at the LHC. Models found by our scan have 
gluinos in the ~ 1.5 — 3 TeV range. This is partially a result of the range scanned, however 
as was scanned up to 2 TeV, the absence of models allowed by our scan with gluino 
masses > 3 TeV, suggests that models with large values of are somewhat more likely 
to run afoul of our constraints, in particular the constraint that the thermal relic density 
be consistent with the WMAP value. 

Gaugino mass unification at the GUT scale is an assumption of this scenario; this 
means that there is a mostly bino LSP with mass in the ~ 250 — 600 GeV range, and a 
mostly wino chargino and neutralino with masses ~ 2mLSP- The Higgsino-like charginos 
and neutralinos generally have masses similar to the wino- like ones. The distribution of 
neutralino masses in these models is shown in Fig. [121 while Fig. 1131 shows the distribution 
of chargino masses; Fig. [14] shows the distribution of the ^ parameter in models found by 
the scan. The distribution of gluino masses is included in Fig. [T5j 

In all models found by this scan, the light sbottom is the lightest squark and is lighter 
than the gluino. The relative lightness of this (mostly right-handed) sbottom (though it 
is not lighter than ~ 1 TeV) is due both to the effect of large Yukawa couplings on the 
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Figure 10: This histogram shows the distribution of BR(b 57) in the models found by 
the scan. The upper and lower limits for BR{b — )■ ,87) which we impose are also shown. 
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Figure 11: This histogram shows the distribution of SUSY contribution to the anomalous 
magnetic moment of the muon (a^ = {g — 2) /2) in models found by the scan. The one- 
loop calculation was done with DarkSUSY; the two-loop contribution was calculated with 
SoftSUSY. 
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Neutralino Masses in Models from the Scan 
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Figure 12: This histogram shows the distribution of neutralino masses in models found by 
the scan. 
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Figure 13: This histogram shows the distribution of chargino masses in models found by 
the scan. 
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Figure 14: This histogram shows the distribution of the Higgsino mass parameter i^t in 
models found by the scan. 
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Figure 15: This histogram shows the distribution of gluino masses in models found by the 
scan. 
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Squark Masses in Models from the Scan 
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Figure 16: This histogram shows the distribution of the first and second generation squark 
masses in models found by the scan. 
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Figure 17: These scatter plots provide a comparison between the bottom squark masses 
and the gluino mass in the models found by the scan. 
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Gluino Mass vs. Light Top Squark Mass 
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Figure 18: These scatter plots provide a comparison between the top squark masses and 
the gluino mass in the models found by the scan. 



sbottom mass RGE and the effect of D-terms, which lower the mass of right-handed down- 
type squarks when M|, > (as we demand in the scan). The lightest stop is also light in 
our models; generally lighter than the gluino. In an appreciable fraction of the models in 
this model set, the heavier stop and sbottom squarks are also lighter than the gluino. First 
and second generation squarks, however, tend to be somewhat heavier than the gluino, as 
is shown in Fig. [TBI 

Since the light sbottom (and generally the light stop and usually the other third- 
generation squarks) are lighter than the gluino, the dominant decay of the gluino will 
be two body decays to either sbottom and bottom quark or stop and top quark. The 
branching fractions for decays to these squarks and quarks are shown in Fig. [191 The 
branching fractions were obtained using SUSY-HIT 1 . 3 |132j . 

The stops and sbottoms produced in gluino decays will in turn decay to charginos and 
neutralinos, including the Higgsino-like charginos and neutralinos due to the large bottom 
and top Yukawas. Then, the charginos and neutralinos will decay either to the LSP and an 
additional Higgs or electroweak gauge boson or to states with taus, as the stau is generally 
the Next to Lightest Supersymmetric Particle (NLSP). Thus the primary signatures of 
the SUSY scenario considered here are cascade decays, with an enriched fraction of b-jets, 
boosted tops, and taus. Clearly this model is discoverable at the 14 TeV LHC. However, at 
the 7 TeV LHC it will be very difficult to discover sparticles in the standard channels. In 
particular, the inclusive gluino pair production cross section in the model with the lightest 
gluino found by the scan (1438 GeV) was investigated at NLO using Prospino 2. 1 |133j . 
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Branching Fraction for Gluino Decays to Third Generation Squarl<s 
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Figure 19: This histogram shows the distribution of branching fractions of the gluinos to 
various third generation squark mass eigenstates and the associated quarks in the models 
found by the scan. 



The gluino pair production cross section for this model is ~ 0.07 fb, while the associated 
squark gluino production cross section is ~ 0.09 fb. This suggests that the production 
cross section for these processes is generically small in this scenario at the 7 TeV LHC. 

6.4 Dark Matter 

In the scenario considered here, the neutralino LSP can have a thermal relic density in 
agreement with WMAP. This is largely due to LSP annihilation through the CP-odd 
Higgs. As can be seen in Fig. [20l ttia is roughly 2mLSP in all of the models selected 
by the scan. As noted above, the lightest stau is often the NLSP; the ratio of the light stau 
mass and the LSP mass is shown for the models found by the scan in Fig. [2TJ We note 
from the distribution of this ratio that there are models found by the scan for which stau 
coannihilation will be important (those where the mass ratio is < 1.1), however in most 
models stau coannihilation becomes negligible. 

We have already noted that efficient annihilation through the CP-odd Higgs and (in 
some cases) coannihilation with the light stau allow the neutralino relic density to be 
consistent with WMAP values. This is different than the case in more standard Yukawa 
Unification scenarios, e.g. [45^150) . where the overproduction of bino-like neutralinos in 
the early universe requires one to introduce a new state, such as an axino, as the true 
LSP. Fig. [22] shows the LSP mass and WIMP-nucleon cross section as compared with the 
XENONIOO bounds. We note that the models found by the scan, which have the lowest 
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CP Odd Higgs Mass vs. LSP Mass 
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Figure 20: This scatter plot shows the distribution of the CP-odd Higgs mass and 
the mass of the neutrahno LSP. That the ratio between the two is generically close to 2 
suggests that neutralino annihilation through the CP-odd Higgs is efficient in all of the 
models selected by the scan. 



Light Tau Slepton Mass vs. LSP Mass 
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Figure 21: This scatter plot shows the the distribution of the lighter f slepton mass and the 
neutralino LSP mass. For models with low values of this quantity {rrif /misp ^ 1-1)) stau 
coannihilation may also play an important role in bringing the thermal LSP relic density 
down to the WMAP value. 
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Figure 22: Here we show the WIMP mass and WIMP-proton scattering cross section for 
models found in the scan as compared to the XENONIOO limits on the spin-independent 
WIMP-nucleon cross section. 

values of spin-dependent WIMP-nucleon cross section, generally have higher LSP masses. 
Conversely models found by the scan with lighter LSPs should be discovered or ruled out 
by the next generation of direct detection experiments. Direct detection constraints were 
important in setting the lower bound on LSP masses found by the scan of ~ 250 GeV. 

7 Conclusions 

We have presented a scenario which adds to the framework of Yukawa-unified SO{10) 
SUSY GUTs a mixing between generations in the charged lepton sector, which produces 
neutrino mixing. We find that the specific values of this mixing angle that one obtains 
from the values of the top, bottom, and tau Yukawa couplings at the GUT scale, in the 
specific 50(10) inspired parameter space considered, suggest sin^ 26*23 ~ 1) in accordance 
with data, in the limit where no other large mixings are present. 

Thus, while this scenario adds a parameter to the framework, the natural value taken 
by this parameter explains observed physics. Another property of this scenario is that we 
find viable points in our parameter space with sparticle masses above the current bounds, 
and quite generally above the reach of the 7 TeV LHC, but likely within the reach of the 
14 TeV LHC. Viable points in our parameter space also have neutralino dark matter in 
agreement with WMAP and direct dark matter detection search constraints. Finally we 
find viable points in our parameter space for which the SM-like Higgs mass is between 
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118 and 122 GeV; a prediction which the 7 TeV LHC will test. In future work, we will 
consider the collider and flavor phenomenology of this scenario in more detail and study 
the importance of the assumption of small neutrino Yukawa couplings made in this work, 
as well as the effects of considering all three generations. We look forward to tests of this 
scenario at the LHC. 
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